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Abstract 

The recent measurement of the muon anomalous magnetic moment shows 
a 2.6cr deviation from the standard model value. We show that it puts strong 
constraints on the parameter space of the two-Higgs doublet model (2HDM) 
II. The dominant contribution of the Higgs bosons comes at the two-loop level, 
and in order to explain the data it favors a pseudoscalar A with a light mass 
range and a large tan/3. At 95% C.L., the upper limit for tua is 29(55) (85) 
GeV for tan/? = 30(45) (60), and tan/? is bounded below at 17. This is in 
sharp contrast to the conclusion one draws from considering one-loop Higgs 
contributions alone. Finally, we also discuss the role of the Higgs contributions 
in the minimal supersymmetric standard model. 
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The recent result on the measurement of muon anomalous magnetic moment (a^) by the 
experiment E821 at Brookhaven National Laboratory has reduced the error to a very 
small level. Comparing the data with the prediction from the Standard Model (SM) 0, one 
gets the deviation 

Aa^ = a;fP - af^ = 426 (165) x IQ-^^ , (1) 

which may suggest the presence of contributions from physics beyond the SM. Taking the 
above numbers at face value, the range of Aa^ at 95% CL. (±1.96(t) is given by 

10.3 X 10"^° < Aa^ < 74.9 x 10"^° . (2) 

Most of the extensions of the SM start with an extended Higgs sector, the simplest of 
which is the two-Higgs-doublet model (2HDM). In particular, the more interesting model II 
0] shares the same Higgs structure as the minimal supersymmetric standard model (MSSM). 
Hence, going beyond the SM to look for extra contributions to a^, the 2HDM should be 
among the first to be examined seriously. 

In this Rapid Communication, we investigate the contributions from the 2HDM II. Strin- 
gent constraints are thus obtained on the parameter space of the model. The results are 
summarized as follows. In the 2HDM II, the dominant contribution actually comes at 
the two-loop level, and in order to explain the data Aa^ is preferred to come from the 
Higgs pseudoscalar A, the mass of which is required to be less than 29 (55) (85) GeV for 
tan/5 = 30 (45) (60) by the 95% C.L. range of Eq.(0). Moreover, tan/? has to be larger than 
17. Recently, there has been some work on the same subject 0]; however, the paper only 
considered the one-loop contributions from the Higgs bosons, which become smaller than 
their corresponding two-loop contributions when the Higgs bosons are heavier than a few 
GeV and substantially smaller for heavier Higgs bosons. It is important to note that the 
two-loop contributions considered here dramatically change the story of the Higgs sector 
contributions to and invalidate most of the results of the one-loop studies. Although 
we focus on the 2HDM II here, a similar conclusion holds for most models with possible 
large contributions from the Higgs sector. In models with flavor- changing Higgs couplings, 
in particular, there are potentially 1-loop contributions substantially larger than the flavor- 
conserving ones 0. However, the Barr-Zee type 2-loop contributions would still have an 
important role to play and should be taken into consideration. This fact has often been 
overlooked in the literature. 

In the MSSM, the dominant contribution comes from the chargino-sneutrino-loop dia- 
grams and it has been shown p that to satisfy Aa^ requires the gaugino mass and the 
smuon mass below about 600 — 800 GeV. Because of additional mass constraints on the 
scalars in MSSM, the total contribution from the Higgs bosons is not at a significant level 
and thus will not affect the conclusion of Ref. . While the typical failure of the MSSM 
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studies to address this kind of 2-loop contribution is a potential problem, our results here 
get rid of the worry, at least for the case of a more generic scalar mass spectrum. 

Many other extensions of the SM have extra contributions to the a^. Some examples 
are additional gauge bosons 0, leptoquarks 0, and muon substructure However, not 
all of them can contribute in the right direction as indicated by the data. Thus, the a^^^ 
measurement can differentiate among various models, and perhaps with other existing data 
can put very strong constraints on the model under consideration. 

Given the mass bound on SM Higgs boson |T^, the Higgs contribution to Aa^ at one- 



loop level is negligible. However, it has been emphasized, in Ref. |Tl| for example, that for 
Higgs boson mass larger than about 3 GeV, the dominant Higgs contribution to actually 
comes from the two-loop Barr-Zee diagram (first discussed by Bjorken and Weinberg) 



with a heavy fermion (/) running in the loop. A m^/m^ factor could easily overcome the 
Q;/47r loop factor. The two-loop scalar contribution with a heavy fermion /, as shown in 
Fig. ID, is given by 

where 

, 1 l-2x(l-x), x{l-x) 

f(z) = -z dx — ^ ^In^ , (4) 

■'^ ' 2 Jo x{l-x)-z z ' ^ ^ 

A^-^ represents the number of color degrees of freedom in /, and Qf its electric charge. Here 
£ denotes a generic lepton, nih is the (scalar) Higgs boson mass, and and A/ represent 
plausible modifications to the Higgs couplings of the fermions (A^ = A/ = 1 in the SM). 
Ref. quoted an electroweak contribution, calculated up to two-loop level, of 

af^ = 152(4) X 10-" . (5) 

Included in this number is a Barr-Zee diagram contribution with a t loop. The numerical 
value of this result barely exceeds the order of 10~^^ (see also Ref. [|13|) and is negative for 



any reasonable value of m^. Moreover, there are other purely bosonic two- loop contributions 



14 1 , in which the SM Higgs boson also plays a role. Nevertheless, all the contributions 



involving the SM Higgs boson are quite small [pt |T3i|T4| . When considering a model with 



an extended Higgs sector, a complete analysis would first require one to subtract the SM 
Higgs contribution and recalculate all the Higgs contributions. This is because the number 
of Higgs bosons, their effective couplings, and mass constraints would be different from the 
SM scenario. Nevertheless, in our study here, we only calculate the Higgs Barr-Zee diagram 
contributions and assume that this does give a very good approximation of from the 
diagrams involving Higgs bosons ( AaJ^'^^^) . Our rationale is as follows. We are interested in 
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the region of parameter space where the Higgs contributions could explain the discrepancy 
of Eq.(|l|), or at least, in the case of models that have other important contributions, do 
play a substantial role in Aa^. Hence, we focus on the region where Aa^^^^^ is at or close 
to the order of 10^^. As we will see below, the possibility of having such largely enhanced 
Higgs contributions comes from the combined effect of coupling enhancements and weakened 
Higgs mass constraints. The coupling enhancement is only to be found among the Yukawa 
couplings, thus illustrating the special importance of the Barr-Zee diagram considered. 

In a model with an extended Higgs sector, we can write the fermion couplings of a neutral 
Higgs mass eigenstate (p^ as 



£f<l>°f 



grnj_ 
2Mw 



2M, 



(6) 



where A/ |^ and Af are the effective scalar and pseudoscalar couplings. The contribu- 
tion of two-loop diagram (Fig. |l|) to is then given by the sum of Eq.(|^) (with rrih = m^o) 
and the pseudoscalar expression 



47r2sin2^^ 



A i 



(7) 



(with m4 = m^o), where 



-z [ dx — r In — ^ 

2 Jo x(l — x) ~ z 



(8) 



Without CP violation in the Higgs sector, Eqs. (|^) or (|^) gives directly the contribution from 
a scalar or a pseudoscalar, respectively. The corresponding contributions with the 7 in Fig. |I| 
replaced by the are suppressed by about two orders of magnitude |Tl[] , hence neglected 
here. We also skip the details about the similar contribution from a charged Higgs boson, 
which involves a boson and is thus suppressed also. Note that a light charged Higgs 



boson less than 80.5 GeV is ruled out by the CERN e^e~ collider (LEP) experiments [[10 
Moreover, analysis of its contribution to 6 ^ S7 leads to a much stronger lower bound — 
380 GeV, as claimed in Ref. [|r^, for instance. 

The 2HDM has three physical neutral Higgs bosons: two scalars h and H, and the 
pseudoscalar A. For model II under consideration, the corresponding nonzero Xf or Af for 
f = t,b, and £(= e, /x, r) are given by 



hiXj) : 

H{Xf) ■■ 
A{Af) : 



cosa 

sin/3 
sina 
sin/3 
cot/3 



sma 
cos/? 
cosa 
cos/? 
tan/3 



sma 
cos/3 
cosa 
cos/3 
tan/? 
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respectively, in the standard notation [Q. What is particularly interesting phenomenologi- 
cally is the enhancement of the couplings of bbH,bbh, HH, i£h,bbA, and HA at large tan/?. 
In fact, the dominant contributions then come from the diagrams with a 6 or r loop. 

It was pointed out in Ref. that the two-loop pseudoscalar contribution to is 
positive while the two-loop scalar contribution is negative in the large tan/3 region. Note 
that this is always true for the dominating contributions with a 6 or r loop, independent of 
the scalar mixing angle a. The reverse happens in the corresponding one-loop contributions, 
but these one-loop contributions are suppressed relative to the two-loop contributions. In 
fact, using the one-loop result in constraining the 2HDM II has been studied extensively 
T6| . Using the two-loop result, however, changes the story dramatically and invalidates 



most of the conclusions from the one-loop studies. 

Applying the Aa^ constraint to the 2HDM II, we need to suppress the scalar-Higgs {h 
and H) contributions, as it comes in the opposite direction as indicated by the data, relative 
to the pseudoscalar contribution. This is in direct contradiction to what is suggested in the 
one-loop studies. At large tan/3 both the dominating contributions from the scalars and the 
pseudoscalar scale roughly as tan/3. For the scalar part, we can adjust the mixing angle a 
to zero such that the contribution from the light Higgs boson h is negligible, and impose a 
large mass hierarchy between the scalar Higgs H and the pseudoscalar A. Then a relatively 
light mass for the pseudoscalar A will give a sufficiently large positive contribution to Aa^, 
or the required Aa^ value could be used to obtain the admissible range for m^. In Fig. |^, we 
show the contribution of Aa^ from the pseudoscalar A versus for various values of tan/3. 
We included the one-loop and two-loop pseudoscalar contributions. The shaded region is 
the 95% CL. range of Eq.(^. The required range of is then given by about 4-29 (15 
- 85) GeV for tan/? = 30 (60). Moreover, a tan/? > 17 is always required. 

What happens when other mixing angles a are chosen? The light Higgs boson h will give 
a negative contribution to and thus offsets the pseudoscalar contribution. Therefore, the 
required range of shifts to a lower value in order to accommodate the data. We show 
the contour plots of Aa^ in the plane of {my^,mh) for a small and a large a at tan/? = 40 
and 60 in Fig. ^ The a = limit, which corresponds to switching off the contribution 
from the scalar h, can be easily read off from the vertical asymptotes. On the other hand, 
a = ±7r/2 corresponds to the maximal contribution from h. The contour plots show how 
the contribution from the scalar h affects the solution to a^. For example, at a = —tt/2 for 
rrih ~ 60 GeV, the 95% CL. required range of lowers to 13-47 GeV for tan/3 = 60. 

At this point, it is interesting to take into consideration other experimental constraints 
on Higgs masses. A collider search for the neutral Higgs bosons in the context of the 
2HDM typically rules out a region of small rrih and m^. In particular, an OPAL analysis 



T7| using the LEP II data up to ^/s = 189 GeV excludes the regions 1 < rrih < 44 GeV 



and 12 < < 56 GeV at 95% C.L., independent of a and tan/3. This is a conservative 
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limit. Details of the exclusion region vary with a and tan/3, and go substantially beyond 



the rectangular box [T^. Essentially, the search for A relies on the process e^e~ Ah, 
therefore, if rrih is so large that this process becomes negligible for A production, there would 
be no limit on m^. In addition, below ^ 5 GeV, the direct search in e~^e^ ^ Ah was not 
included because the detection efficiency vanishes and the total width only provides very 
limited exclusion. The OPAL exclusion region is roughly sitting at the center on the axis 
going up to about m/j = 60 GeV in the plots of Fig. |^, and cuts out part of the admissible 
region of the Aa^ solution. For rrih larger than the OPAL limit the admissible range is 
roughly from 3 to 50 GeV for tan/5 = 40. At larger tan/? (60 as the illustrated example), 
the range widens, especially at the upper end; but a middle range (about 4 to 15 here) is 
lost as Aa^ gets too large. There is another admissible small window at very small rrih (a 
few GeV) with large m4 . This region is indeed dominated by the one- loop contribution from 
h, but even here, the two-loop contribution (from h) has an important role to play. This 
tiny ruh window and the similar solution with this kind of small rriA are in fact excluded by 



Upsilon decay and some other processes [Q. 

In addition to the constraint from direct search, there are also other constraints on the 
masses of the Higgs bosons coming from the electroweak precision data. While a comprehen- 
sive treatment of the topic is really beyond the scope of the present study, we discuss below 



the basic features, using results from a recent paper [T^. The scope of the latter study is 
limited to the large tan/? region and a = /5 ~ 7r/2. There, the ratio mh/rriA is constrained 
at 95% C.L. (based on the Bayesian approach) via the function 



<] >_ ( 39 y 

\tan(3 J 



ml 



G 
where 

Solving for mh/rriA at tan/? = 60,45,30, we obtain, respectively, 

0.3 < < 3.2 

— niA — 

0.2 < ^ < 5.1 

— niA — 

0.07 < < 14.4 . 

— niA — 

We can see that the precision electroweak data prefers a region close to the diagonal of the 
m/i verses plot. If one naively imposes the result of 0.3 < ^ < 3.2 at tan|!/5 = 60 
onto Fig. 3(d), which is at a different a but with which a similar result is expected to be 
valid, together with the direct search limit nih ^ 60 GeV and the requirement, only a 
small "triangle" is left. This triangle is bounded by ruh = 60 GeV, mh/mA = 3.2 and the 
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contour of = 10 x 10"^'^ [labeled by "10" in Fig. 3(d)]. The surviving parameter space 
region, however, is in the more favorable "larger" mass area. In particular, it re-enforces 
our previous comment at the end of the last paragraph that the one-loop dominating tiny 
rrifi window of solution to is ruled out. The situation for the other a values is expected 
to be similar. It will be very interesting to have the complete phenomenological analysis 
combining all the constraints on the 2HDM. 

Finally, we discuss the role of the Higgs sector contributions to in the MSSM. The 
LEP bound on the Higgs boson masses is in the range 85 - 95 GeV ||20|. From the above 
result, one may naively conclude that if a Higgs particle is just around the corner, it could 
have an important role to play in a^. However, there are some strong theoretical constraints 
on the relation of the Higgs boson masses in the MSSM. At the large tan/3 value required, 
one has < ~ Most of the Higgs contributions to cancel among them- 

selves. Moreover, a small Higgs boson mass may require a /i parameter so small that the 
chargino/neutralino contributions to get far too large. In fact, we have checked and found 
no interesting solution within MSSM in which the Higgs contributions play a substantial 
role. In the admissible range of chargino/smuon masses found in Ref. 0, the total Higgs 
contribution is only about 1% of the SUSY contribution. The importance of this null result 
should not be underestimated. 

We conclude that the new measurement on the muon anomalous magnetic moment con- 
strains severely on the parameter space of the 2HDM II, and our results, including one-loop 
and two-loop contributions, change dramatically from the conclusion that one draws by 
using only one-loop results. 
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Figure Captions :- 

Fig. 1 — The dominant two-loop graph involving a scalar or a pseudoscalar boson that 
contributes to a^. 

Fig. 2 — Aa^ (one-loop and two-loop summed) versus the pseudoscalar Higgs boson mass 
ttia for various values of tan/5. 

Fig. 3 — Contours of Aa^'^^^ in unit of 10~^° in the plane of (m^, mh) for (a) a — —n/S and 
tan/3 = 40, (b) a — — 7r/2 and tan/3 = 40, (c) a — —n/S and tan/3 = 60, and (d) a — — 7r/2 
and tan/3 = 60. 
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FIGURES 




FIG. 1. The dominant two-loop graph involving a scalar or a pseudoscalar boson that con- 
tributes to a^. 
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FIG. 2. Aa^ (one-loop and two-loop summed) versus the pseudoscalar Higgs boson mass rriA 
for various values of tan/3. 
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FIG. 3. Contours of Aa^'^^*^ in unit of 10^^'' in the plane of (m4,m/j) for (a) a = — vr/S and 
tan/3 = 40, (b) a = — 7r/2 and tan/3 = 40, (c) a = —vr/S and tan/3 = 60, and (d) a = — 7r/2 and 
tan/3 = 60. 
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